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DIRECT INTEGRATION TRANSMITTANCE MODEL 

V. G. Kunde 
W. C. Maguire 

Laboratory for Planetary Atmospheres 
ABSTRACT 

A transmittance model has been developed for the 200-2000 cm ^ region 
for interpretation of high spectral resolution measurements of laboratory 
absorption and of planetary thermal emission. The high spectral resolution 
requires transmittances to be computed monochromatically by summing the 
contribution of individual molecular absorption lines. A magnetic tape atlas of 
H2O, O3, and C 0 2 molecular line parameters serves as input to the trans- 
mittance model with simple empirical representations used for continuum 
regions wherever suitable laboratory data exist. The theoretical formulation 
of the transmittance model and the computational procedures used for the 
evaluation of the transmittances are discussed, and application of the model to 
several homogenous path laboratory absorption examples is demonstrated. 
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1. INTRODUCTION 


Improved spectral resolution in infrared measurements has required a cor- 
responding improvement in the theoretical and analytical techniques for interpre- 
ting the measurements, particularly in the area of molecular absorption. Data 
of intermediate to high spectral resolution in the infrared are now available for 
the Earth from the Nimbus 4 interferometer spectrometer (HANEL et al. ^); 
for the Moon, Venus, and Mars from ground-based observations (CONNES et 
al. HANEL et al. ( 3 ), and RANK W); and for Mars from the interferometer 
spectrometer on the Mariner 9 orbiter (HANEL et al. ( 3 )). High spectral res- 
solution infrared measurements are also available for stellar sources (CONNES 
et al/ 6 ), GEBALLE et al/ 7 )). Representation of the molecular absorption by 
band models and empirical fits to laboratory absorption data is no longer 
adequate for interpretation of the high spectral resolution measurements 
(DRAYSON ( 8 )). The need for high spectral resolution and for versatility in 
changing parameters such as spectral resolution, instrument function, and 
molecular line parameters requires new techniques for handling molecular 
absorption. 

One of the techniques that has been successfully developed is the direct inte- 
gration method, which involves the computation of the monochromatic absorp- 
tion spectrum by numerically summing the contributions of the individual 
molecular absorption lines. This technique has been applied to interpretation 
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of planetary and stellar infrared spectra (DRAYSON ( 8 ), KYLE ( 9 ), KUNDE ( 10 ), 

AUMAN ( 11 ), QUERCI et al. ( 1 ^)). Direct integration programs for homogenous 

path applications have been developed by CALFEE and by ARNOLD et al. ( 1 ‘ 1 ). 
This paper describes a direct integration transmittance model developed by the 

authors for interpretation of laboratory absorption and planetary thermal 

emission spectra in the 200-2000 cm - * region. The application of the model to 

laboratory absorption data is also presented. The model has been used previously 

to interpret planetary thermal emission spectra of the Earth and Mars obtained 

with the interferometer spectrometers on Nimbus 4 and Mariner 9 (Hanel, et al/ 5 \ 

Kunde, et al/ 15 )). In subsequent papers the transmittance model will also be 

applied to the interpretation of ground-based planetary observations of Venus 

and Mars . 

The formulation of the infrared radiative transfer equation for homogenous 
and slant path problems is given in Section 2; the formulation of monochromatic 
molecular absorption coefficients is outlined in Section 3; and the sources of 
the various molecular line parameters are discussed in Section 4. 

2. INFRARED RADIATIVE TRANSFER THEORY 

This section reviews the general theory and mathematical formulation 
governing the radiative transfer of infrared radiation through a single layer of 
absorbing gas in which the physical state of the gas does not vary (homogenous 
path problem), and a series of layers at varying temperatures and pressures 
(slant path problem). The homogenous path problem is generally applicable to 
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laboratory absorption measurements, while the slant path problem relates to the 
transfer of radiation through a planetary atmosphere. A specific application of 
the formalism developed in this section and the adopted numerical techniques 
will be discussed in Section 3. 

For computational efficiency, the general procedure adopted for both prob- 
lems is to compute the average transmittance of radiance for Av intervals which 
are much smaller than the instrumental resolution. The observed transmit- 
tances or radiances are then obtained by convolving the computed spectrum with 
the instrument function. In this paper, the value Ap = 0. 1 cm -1 has been 
adopted. 

Homogenous Path Theory 

With the assumption of local thermodynamic equilibrium, the spectral radi- 
ance (N„) from a slab of a mixture of absorbing gases at uniform temperature 
(T) and pressure <P) may be written in its most general from as 

N„ = N®ex p [-S k >.T)u.] + B (> (T)[l-exp(-Xk>.T)u i ) ] (X) 

i * 

with the optical path length u = pL. The geometrical path length is L, the total 
molecular absorption coefficient for the i-th gas is kl, (P, T), the gas density 
is Pj , the Planck function is (T), and N° is the external incident intensity on 
tiie slab of gases at L = o. The first term on the right represents the transmitted 
radiance; and the second term, the thermal emission of the gas. If the thermal 
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emission is negligible (ku very small), Equation (1) reduces to 


r ,=-> p u i] 


where t v is the monochromatic transmission through the absorbing gas. The 
transmittance averaged over a small spectral interval, Av, is 



which is numerically evaluated using Legendre -Gauss quadrature (STROUD and 
SECREST (16)). Each Ai> spectral interval is divided into a variable number of 
subintervals (k = 1, 2. . . K) with four-point Lengendre-Gauss quadrature ( C = 


1, 2. . .4) used for each subinterval. Equation (3) now may be written as 



( 4 ) 


where 



( 5 ) 


with and v A being the upper and lower wave number boundaries of the sub- 
interval, respective; and wg and xg the Legend re -Gauss weights and abscissa, 
respectively. 
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Associated with all spectrometers is an instrument function, f( |r - y n I, a), 


which degrades the true monochromatic spectrum at wave number v n to 


f(l*'-»' n l> a ) Vb 


where a represents the full width of the instrument function at half maximum. 
The instrument function has been normalized 



+ oo 

f(b-»' n l> a )d>'= 1. 


(7) 


Assuming Ay is small enough that f( \ v ~ % I, a) is essentially constant over 
Ay , Equation (6) can be approximated by 

V£ 7(l ^'’ a > 7 - A ' (8) 

with f( \v — v n \, a) representing the average value of the instrument function over 
the A v interval. The principal advantage of evaluating the instrument -averaged 
transmittance from the approximation in Equation (8), using the transmittance s 
initially averaged over Ay intervals rather than directly from Equation (6), is a 
large saving in computer storage. In the approximate scheme the monochro- 
matic transmittances need not be save after t^ v is computed, thus saving a 
factor of ~100 in core storage. The approximation scheme is adequate as long 
as a > Ay, 
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The convolution of Equation (8) can also be represented in terms of the 
Fourier transforms of the functions involved: 

r*(t) =7*(t) -7* p (t) (9) 

with t indicating the time domain and the asterisk denoting (he Fourier trans- 
form. In the time domain the convolution reduces to a multiplication. For a 
large spectral range or for small values of hv it becomes more efficient to de- 
termine r„ through the fast Fourier transfrom technique rather than through 
the direct convolution process. 


Slant Path Theory 

The emergent thermal emission spectrum at the top of the atmosphere has 
its origin in the radiant energy processes occurring in the planetary atmosphere 
and at the planetary surface. The emergent monochromatic radiance, neglect- 
ing scattering processes, is given by the radiative transfer equation as 

N ,( S T) = %B i ,(T 8 )r (1 (0,S T ) + r F- r v (0,S T ) + B„[T(s')] dr/s', S T ). (10) 

•'t (o,o) V ; 


The integration is through the planetary atmosphere along any geometric slant 
path (s) from the planetary lower boundary surface (s = 0) to the top of the 
atmosphere (s = S T ). A schematic of the radiation components and definition of 
some of the slant path parameters is shown in Figure 1. The emergent radiance 
consists of the surface emission, the downward atmospheric radiation flux re- 
flected from the planetary surface, and the atmospheric emission. The emissive 
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characteristics of the surface are described by e v atmospheric transmittance 
between levels s and S x by r v (s, S T ), reflectivity for a perfectly diffuse re- 
flector by r„ , and the downward Atmospheric flux by 7rF~ . The subscript s 
denotes planetary surface parameters. The reflected component can generally 
be neglected as most surfaces are blackbody or near-black-body in character,. 
(r„ ~ 0). An exception to this occurs over surface areas where the strong 
reststrahlen bands may yield a significant reflected component; 

i t ^ - 

The monochromatic transmittance through an inhomogenous medium of a 
mixture of gases is 

S T ) = exp ^ k v ' (P, T) p.ds! J. (11) 

i 

i 

With the optical path length for the i-th gas 



( 12 ) 


we can approximate the integration along the atmospheric slant path by a series 
of isothermal, constant pressure layers (j = 1, 2 , . . . , J) as ; 


r y (s, S T ) = exp C (Pj- ty ^ ] 

) i 


(13) 


where the j -summation includes all atmospheric layers between levels s and 
S T , and P and T are the average pressure and temperature for a layer. We 
can define a slant path average transmittance, similar to the homogenous path 
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average transmittance, in terms of fou'r^point Legendre^auss quadrature yioKi;, 



The radiative transfer formulation given above describes only the general >■. 

approach used to determine the emergent spectrum for any atmospheric slant . 'j 
path, s, through a model atmosphere?; A mote detailed description of. the. ' i 
techniques used to relate the slant path, s, to the local vertical path, h, for a 

• ■'.?*' ■ * , ( t. ’ . *. ' ■ ; ' iy't 'ri i* 

spherical atmosphere and for the determination of the optical path lengths of 
the various absorbing gases along s has been given previously (KUNDE ( 18 )). 

The emergent spectral radiance in a finite spectral interval, Av , is 


n a„( s t )-/* N,(S x )d* 


'* ' t". : 


i 


j* i* 


(15) 


With Equation (10) and assuming that Av is small enough to allow the spectral 
quantities to be evaluated at the midpoint of Av (denoted by v'). Equation (15) 
may be expressed in terms of the average transmittance as 




Ap 




T i^ 0, S t^ 


; V W t a,(°> S T ) + V *V S T ) Vi B y ,(T(s'))dF Au (s', S T ). (16) 

' * " Av(0, 0) 


We can numerically approximate Equatioii (16) by 


N a„(S t ) 


* »r ■* '• til .■<< 


av v t / _ _ r — — i 

"2P = &r) + s t) + [b.'W 5 ')) s T ) J = i . .a i tm. v. 
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where the atmosphere has been divided into a series of isothermal layers with 
the integral term of Equation (16) being replaced by a sum over the contributions 
of all atmospheric layers. The difference in transmission across the layer is 
A^Ap- w ith the layer increment taken as As' = 0. 1 km. To evaluate Equation 
(17) the average transmittance 7 (s, S T ) is precomputed for a set of atmos- 
pheric levels (s p , p = 1,2,. . . N), with N being approximately 30. The N levels 
are selected to yield approximately equal increments of A r. The vertical layer 
mesh usually varies from 0. 5 to 4 km with the smaller layers being concentrated 
near the lower surface. The emergent spectrum for each Av = 0. 1 cm' 1 interval 
is computed from Equation (17) for As' = 0. 1 km steps with the average trans- 
mittance obtained by interpolation from the precomputed transmittance array, 

7a„ (Sp, S T ). Computing the average transmittance on a coarse vertical scale 
is essential for making the direct integration calculations tractable. 

In order to compare the theoretical outgoing radiances with measured radi- 
ances (N ), we must convolve the theoretical radiances with the instrument 
n 

function, f( \v - v n |, a), 

(18) 

This convolution is accomplished by a fast Fourier transform procedure in the 
same manner as described for the homogenous path convolution. 
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The numerical algorithm for determining the emergent radiances incorpor- 
ates the following steps: 

(1) Specify atmospheric model, vertical temperature profile, vertical gas 
concentration, and surface pressure 

(2) Specify molecular absorbing gases and appropriate line parameters 

(3) Set up vertical matrix of atmospheric temperatures, atmospheric 
effective pressures, and optical path lengths for each absorbing gas 
for each atmospheric layer 

(4) Determine the average transmittance for a 0. 1 cm' 1 interval for J 
atmospheric levels from Equations (13) and (14) 

(5) Determine the emergent radiance for each 0. 1 cm’ 1 interval from 
Equation (17) across the spectral region of interest 

(6) Determine the emergent radiance averaged over the instrument function 
as specified in Equation (18), for comparison to the observed spectrum. 

This algorithm has been programmed entirely in FORTRAN IV for an IBM 360 
operating system. 

3. MONOCHROMATIC MOLECULAR ABSORPTION COEFFICIENT 

FORMULATION 

In this section we will describe the determination of the monochromatic 
molecular absorption coefficient at the wave number mesh points necessary 
for the monochromatic transmittances. 
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A molecular absorption line can be described by four pai'ameters: line 
position (v Q ), integrated line intensity (S), relative line shape (b), and line half- 
width (a). The absorption coefficient for a single line is 

% = S b(y - v Q> ft) (19) 

where the relative line shape b(t> - v 0 , a) depends on the type of line-broadening 
mechanism considered. 


The two line shapes of importance in the terrestrial atmosphere are the 
thermal Doppler 

R ^ d ) 2 2"2 


V£n2 

b (,- V « D )=— — ex p 


Of D n 


( 20 ) 


and the collision-broadened Lorentz 


“c 1 

b(i> - v a ) = — — — ~ 

Jr (y - v G y + a c 2 


( 21 ) 


with a D and a Q being the Doppler and collisional line half-widths, respectively. 
The mixed Lorentz -Doppler line shape is 


/ 2n 2 1 / x 

b(a,u ) = / H(a, u) 

M n a D 


( 22 ) 


where a = (a c /a D ) y/Jrii and u = [ (v - r 0 )/« D ] ,/En 2. The function H(a, u) has 
been evaluated following a numerical scheme of YOUNG Kinetic theory 

predicts the temperature and pressure dependence of the collisional half -width 
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to be 


a c (P,T) = a c (P r ,T f )-(-^) (23) 

with 7 = 0. 5. The subscripts r refer to reference conditions. For a mixture of 
gases, the pressure P is usually converted to an effective pressure, ^P e ^, 

P e = P + (B-l) Pa (24) 

to correct for the influence of self -broadening on the line half-width. The partial 
pressure of the absorbing gas is P a and B is the self -broadening coefficient. 

The collisional line half -width now may be written as 

a c (P, t TI = «c(P [ , T 1 )^)- (25) 


The molecular absorption coefficient at a specific wave number is the sum- 
mation of the contribution from each of the individual molecular absorption lines, 


(P,T).y's„,(T)b(.-. 0 ,«„<P.T)) 

m m 


For planetary atmospheres the monochromatic absorption coefficient varies by 
orders of magnitude on a scale of 1/1000 cm -1 for wave numbers near the line 
center, while in the line wings the absorption coefficient changes much more 
gradually. To retain a high degree of accuracy in the transmittance computa- 
tions, the rapidly varying behavior of the absorption coefficient requires a wave 
number mesh spacing of the order 1/1000 cm -1 near the line center, while 
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1/10 cm -1 will suffice in the line wings. The method used to evaluate Equation 
(26) is patterned after the numerical technique developed by DRA.YSON 
Basic to this procedure are a division of each line into direct and wing contri- 
butions and the establishment of a variable wave number mesh for the evaluation 
of k v (P, T). 


The direct contribution includes all wave numbers within a distance 5 of the 
line center, and the wing contribution includes all wave numbers equal to or 
greater than 6, 


D W 

« = £ + £ 

V &v<8 ii/>8 


(27) 


where 8 has been chosen as 3. 5 cm 1 . The total molecular absorption coeffi- 
cient now becomes 


m m 


(28) 


which can be rewritten 


d w 

k s ,(P e ,T) = k,(P e ,T) + k ! , (P e ,T) 


with 


(29) 


C< P e> T )=53 fi I( P e» T ) I*"' 

m 

k >e- T ) = S»> e .T). IAH<5 


(30) 
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The technique of separating each line into direct and wing contributions has the 
advantage that one need consider only lines in the interval v - 5 to v+ 5, rather 
than all lines, for computing the direct contribution, thus insuring a considerable 
reduction in the amount of computer time required. A second advantage is that 
the wing contribution may be easily adjusted for empirical fitting to observed 
spectra in spectral regions where the theoretical treatment is not adequate. An 
example of this will be discussed later for the water vapor continuum in the 500 
and 1000 cm -1 windows. 

The determination of the wave number mesh for the direct contribution will 
be described briefly with the aid of Figure 2. Cases where no line, one line, 
and two lines fall within the Ai> = 0. 1 cm -1 interval are shown in Figure 2(a), 

2(b), and 2(c), respectively. The circled numbers indicate the subinterval 
number; the dashed lines, the subinterval boundaries; and the short lines near 
the bottom of each figure indicate schematically the wave number mesh. The 
wave number scale is distorted near the line centers for clarity of presentation. 
The distance d, from the line center is subdivided into subintervals, each con- 
taining four mesh points whose positions within the subinterval are specified 
by the Legendre-Gauss quadrature abscissa points. The distance dj is usually 
chosen as 0. 01 cm -1 and the number of subintervals in dj is a variable in the 
range of from one to five with the higher number of subintervals yielding more 
accurate transmittance values. The number of subintervals is usually taken as 
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two, which allows the transmittances to be computed with a numerical accuracy 
of two to three significant figures, an accuracy which is sufficient for most 
applications. After the subintervals are established near the line center, the 
remaining sections of the interval, d 2 and d 3 in Figure 2(b), d 2 » d 3 , and d 4 in 
Figure 2(c), are each assigned as one subinterval. For more than two lines in 
a 0. 1 cm -1 interval the division scheme of Figure 2(c) is applied repeatedly. 
Without lines in the 0.1 cm" 1 interval, the four-point Legend re -Gauss quad- 
rature is applied to the entire 0. 1 cm" 1 interval. The scale of the variable 
wave number mesh is thus of the order 0. 001-0. 01 cm -1 . 

The mixed Doppler- Lorentz line shape is used in the direct contribution if 
the total pressure is less than 100 mb and the distance from the line center is 
less than 0. 2 cm' 1 . The Lorentz line shape is used for the remaining direct 
contribution. 


The collisional line shape is applicable for the wing contribution which, 
from Equations (19), (21), and (30) and using the condition A v >a c (P r , T r ) 


w o 

k (P , T) = k (T)- 


(31) 


with 


VE 


(r.T,) /TX 7 


m * (y-v 0 y 


(32) 
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The function k° (T) changes slowly with respect to wave number and thus may 
be precomputed at a coarser wave number grid 0.5 cm -1 ) than required for 
the direct contribution with intermediate values obtained by interpolation. 


From Equations (19), (29), and (31), the total monochromatic absorption 
coefficient necessary for determining the monochromatic transmittance for the 
i-th gas and j-th layer (see Equation 13) is 


K.1 m Til J 



e. 

J 


(33) 


The molecular line parameters for the evaluation of Equation (33) are stored on 
magnetic tape, ordered by ascending wave number, with the parameters v Q , 
Icode, « c , 7, S(175 K), S(200 K), . . . , S(300 K) preserved for each line and 
with Icode identifying the species of gas. The present atlas contains ■"-6000 
H 2 0 lines, ~5000 0 3 lines, and -TO, 000 lines of C0 2 in the 200 to 2000 cm* 1 
range. The sources of these line parameters and the modifications for em- 
pirical line shapes are the subject of the next section. 


4. SOURCES FOR MOLECULAR TRANSMITTANCE PARAMETERS 
The molecular line parameters for carbon dioxide, water vapor, and ozone 
were collected from a variety of sources. The molecular parameters for C0 2 
are generated in the following manner. The line positions and strengths have 
been determined for 68 bands of C 1 2 0 2 16 , using standard spectroscopic formu- 
lation for a linear molecule, from band center positions, energy levels, 
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rotational constants, and band strengths provided by BENEDICT (20). The only 
exception is that the band intensities of DRAY SON (“^ are used if significant 
differences exist. The band parameters for the stronger C 12 0 2 16 bands are 
listed in Table 1. Also included in the computations are 70 isotopic bands of 
carbon dioxide. The strengths of these bands relative to the corresponding 
bands of C 12 0 2 16 are taken to be the same as the abundance of the species 
relative to C 12 O, 16 . The relative abundances are calculated from the ter- 
restrial abundances of the isotopes of carbon and oxygen. 

The self- and nitrogen-broadened line half -widths for carbon dioxide as a 
function of rotational quantum number were obtained from YAMAMOTO et al. ( 22 ), 
and were assumed to be valid for all vibrational bands. 

The C0 2 continuum is substantially sub-Lorentzian in the 1000 cm -1 re- 
gion. BIGNELL et al. ( 23 ) have reported experimental work from which they 
concluded the absorption by C0 2 outside the range from 560 to 790 cm" 1 is 
negligible. WINTERS ( 24 ), following suggestions by BENEDICT, found that an 
exponential modification of the Lorentz line shape beyond some minimum dis- 
tance from the line center improved the agreement of calculated and experi- 
mental absorptance for CO, in the 2400 cm" 1 region. Recently, BIRNBAUM ( 2 ^) 
has developed a theoretical basis for a line shape which is Lorentzian near the 
line center and exponential sufficiently far from the line center. The point at 
which the shape becomes exponential depends on the duration of collision. 
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BURCH has found that the wings of nitrogen-broadened CO, lines absorb 
only about 1 percent as much in the 780-900 cm" 1 regions as if they had the 
Lorentz shape. To fit the experimental C0 2 continuum we have assumed an 
exponential modification of the Lorentz line shape 

b = b Lorentz t 1 ' ' *0 1 ^ ‘'min 

r bi (34) 

b = b Lorentz eX P ["* ( |t '" ^ “‘'min) J I* " *0 1 > "min 

beyond v min =3.5 cm -1 from the line center for all C0 2 lines in the 500 to 800 
cm -1 region. The parameters a and b have been calculated using WINTERS' 
et al. values a = 0. 08 and b = 0. 8 as starting values and running a mesh of 
a's and b’s to find the parameters which reproduce the experimental continuum 
absorption. The best fit is obtained with a = 1. 4 and b = 0. 25. The above ex- 
ponential correction was also used for the C0 2 continuum in the 500 cm -1 window 
region. A fit to the self-broadened C0 2 continuum of BURCH < 26 ) yields values 
of a and b approximately the same as the nitrogen-broadened continuum values. 

The accuracy of the theoretical calculations was checked by computing the 
transmission for various homogenous path pressures and optical masses and 
comparing with published experimental spectra (BURCH et al. ( 27 )). Two of 
these cases are shown in Figure 3, with the theoretical transmission being 
computed with a triangular instrument function of 4 cm" 1 full width at half- 
maximum. The solid line represents the measurement, and the dashed line, 
the calculation. The agreement is better than 10 percent for both cases. The 
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errors in the strong Q-branch regions are enhanced by imprecise knowledge of 
the true instrument function and do not reflect a major discrepancy between the 
theoretical and experimental transmittance s. However, the theoretical com- 
putations need improvement for the pressure-broadened line shape within 5 cm -1 
from the line center. The investigations of BURCH et al. have indicated 
significant deviations from the Lorentz shape for self-broadened C0 2 lines in 
the 2400, 3800, and 7000 cm' 1 regions. Experimental investigations are needed 
to determine the deviation for N 2 -C0 2 and C0 2 -C0 2 collisions for the 667 cm -1 
band. 

The line data used for water vapor are taken from BENEDICT and 
KAPLAN < 29 > for the pure rotational lines and from BENEDICT and CALFEE < 30 ) 
for the 1595 cm' 1 band. 

One of the most difficult problems to treat in atmospheric transmission 
models is the water vapor continuum regions. The difficulties arise because of 
uncertainties in the effect of aerosols, in knowledge of molecular line shapes in 
the far wings of lines, and in knowledge of the physical mechanisms responsible 
for the water vapor absorption. Water vapor continuum absorption in the 500 
cm" 1 and 1000 cm' 1 windows was discussed in the Nimbus 3 IRIS paper 
(CONRATH et al. The water vapor continuum absorption coefficients 

adopted for the Nimbus 3 investigation gave good agreement in the 1000 cm' 1 
window, but overestimated the absorption in the 500 cm -1 window. 
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Since the publication of Reference 31, a number of pertinent laboratory and 
free atmosphere investigations concerning water vapor continuum absorption 
have been published. Laboratory measurements by BIGNELL ( 32 ) and 
BURCH indicate the water vapor continuum absorption in the 500 cm" 1 and 
1000 cm" 1 regions consists of two components, one proportional to the total 
pressure and the other to the water-vapor partial pressure: 

P'h^o’ T) - kj (T) P + k 2 (T) p H ^ 0 (35) 

BIGNELL ( 32 ) has also reinterpreted, in a semiquantitative sense, previous 
atmospheric and laboratory measurements in these spectral regions and con- 
cludes that the previous work in general substantiates the two -component 
mechanism. Recent interpretations of atmospheric measurements by 
PLATT, < 33 > 34 ) WARK, ( 35 > and HOUGHTON and LEE (36) have further veri- 
fied the Bignell and Burch results. 

The two-component continuum has been incorporated in the transmittance 
model; the adopted kj and k 2 values are illustrated in Figure 4. Laboratory 
values of k j in the 1000 cm -1 regions are not available due to experimental 
difficulties. The laboratory values of Burch for this region shown in Figure 4 
represent only the upper limits to k j. The investigations of WARK ( 33 ) and 
HOUGHTON and LEE ^ 33 ^ gave satisfactory results with the k t value in the 
1000 cm' 1 window essentially zero. Based on their success, we have forced 
k, to approximately zero by extrapolation of the measured k, component from 
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the 500 cm -1 window region and of the computed Lorentz foreign-broadened 
wing component of the 1595 cm' 1 H 2 0 band. The k t and k 2 components at 175 K 
were scaled from the 300 K values using the temperature dependence of Bignell. 

Also shown in Figure 4 is the wing contribution of the rotational water-vapor 
lines using the modified Van Vleck-Weisskopf (M-VVW) line shape described 
by FARMER (37). All self -broadened components have been computed for a 
self -broadening coefficient of 6. 3. As illustrated in Reference 31, no theoret- 
ical line shape is adequate for the 500 cm -1 and 1000 cm -1 window regions. 
However, M-VVW line shape gives fair agreement with the measurements of 
PALMER ( 38 > in the 250-350 cm -1 range and, for continuity, the adopted k t and 
k 2 components are normalized to the M-VVW at 400 cm -1 . The Lorentz line 
shape is used for the 1595 cm' 1 H 2 0 band with the normalization to kj and k 2 oc- 
curring at 1400 cm' 1 . The normalization results in slight discontinuities, 
which, however, create no problems in practice. The water vapor continuum 
absorption in the 400 to 1400 cm' 1 range thus is based on laboratory data com- 
puted from Equation (35) with the kj and k 2 values determined numerically. 

For water vapor in this region, Equation (35) replaces the term k^ (P e , T) of 
Equation (29). 

The adopted water vapor continuum absorption coefficients are employed 
in the transmittance model in an empirical sense, to be used as part of a work- 
ing model until more knowledge is acquired of the opacity sources and their 
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governing physical mechanisms. While the empirical model may still need 
considerable adjustment with regard to the actual numerical magnitudes of k, 
and k 2 , it is felt that the two-component concept, especially the p H 0 term, is 
an important step forward in coping with water vapor continuum transmittance 
problems. 


Homogenous path comparisons for the 1595 cm -1 H 2 0 band are shown in 
Figure 5. The experimental data of BURCH et al. varied in spectral resolu- 
tion from ~5 cm -1 at 1200 cm -1 to ~20 cm' 1 at 2000 cm -1 . Thus the R branch 
side of the band does not exhibit much line structure. The agreement in the P 
and R branches is better than 5 percent, with the theoretical transmittances 
being slightly too transparent in the continuum. 


The first-generation ozone line parameters was not very satisfactory, as 
evidenced in the radiance comparisons of Reference 31. A second-generation 
set of ozone lines is now available; the main improvement is the inclusion of 
numerous upper-state bands. This line set seems to be more satisfactory 
based on comparisons to balloon flight measurements by GOLDMAN et al. ( 39 ) 
and to laboratory data by AIDA (40 >, Included in this line set is the v 2 band at 
701. 4 cm 1 , a band which is of importance to atmospheric temperature inver- 
sions in the 667 cm -1 C0 2 band region. Theoretical calculations for this band 
are compared to the laboratory spectra of McCAA and SHAW < 41 ) in Figure 6. 
The number of lines considered for the v 2 band was 2158, with a collisions! 
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half -width of 0. 07 cm -1 at 296 K and a Lorentz line-shape, Self-broadening 
effects were assumed to be unimportant. The general agreement is approxi- 
mately 5 to 10 percent which is satisfactory for initial applications to the tem- 
perature inversion problem. Absorption by these lines is demonstrated for a 
slant path in the paper on application to the Nimbus 4 IRIS spectra (KUNDE ( 15 )). 

5. CONCLUDING REMARKS 

A direct integration transmittance model has been developed for 
interpretation of spectra of intermediate to high spectral resolution recorded 
in the laboratory or in the atmosphere. This procedure, though time-consuming 
on present digital computers, avoids limitations and approximations necessary 
when employing empirical representations and band models for computation of 
molecular absorption. The main advantages of the direct integration technique 
are the ability to attain high spectral resolution and to easily vary this resolu- 
tion; to facilitate different types of line shapes readily, including mixed line 
shapes; to handle line wing effects; to cope with temperature, pressure, and gas 
concentration variations, and to avoid assumptions necessary to reduce a slant 
path absorption to the equivalent homogenous path absorption. The central con- 
cepts for making the direct integration a tractable problem are the variable wave 
number mesh and the division of each line into direct and wing contributions. As 
the numerical procedures are essentially exact, the accuracy of the predicted 
transmittances is limited by the accuracy of the input molecular parameters. 
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As a consistency check on the input line parameters and the transmittance 
model, theoretical predictions have been compared to existing laboratory data 
for several representative homogenous path cases. The comparisons gave 
transmittance disagreements of 5 to 10 percent, a number which is probably in- 
dicative of the absolute accuracy of present state-of-the-art molecular 
transmittance s. 

The accuracy of the slant path model and the predicted transmittances has 
been assessed by comparison to Nimbus 4 IRIS data, and Mariner 9 IRIS data, 
with both applications yielding agreements of 5 to 10 percent in emergent 
radiances through the 200 to 2000 cm -1 region. The model is presently being 
tested further against high spectral resolution terrestrial transmittance spectra 
( ~ 0. 2 cm -1 ) and downward sky emission spectra in the 450 to 1400 cm' 1 re- 
gion obtained from McDonald and Mauna Kea Observatory. 
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Table 1 


Line Parameters and Band Intensities for the Strongest C 12 0 2 16 Bands at 300 K. 


Band No. 

Level 

Lower State 
Energy 
(cm"* 1 ) 

Band Center 
(cm -1 ) 

Band Intensity* 

Lower 

Upper 

cm -1 (atm cm) -1 300 K 

cm" 1 (mol cm -2 ) 

1 

(10°0,02°0) 1000 

0301{ll 1 0,03 1 0) n 

1388.187 

544.283 

0. 10 b 

2.4 x 10” 22 

2 

(ll 3 0,03 3 0) n 0303 

0402(12 2 0,04 2 0) n 

2003.309 

581.697 

0.0084 b 

1.5 xlO -22 

3 

0202 

OSOlflA, 03 1 0) n 

1335.129 

597.341 

0.14 

3.4 x 10" 21 

4 

0101 

0200(10 0 0,02°0)^ 




-19 




667.379 

618. 033 

4.27 

1.43 x 10 

5 

(10°0,02°0) n 0200 

OSOlfllVoa^jj 

1285.412 

647. 058 

0. 591 b 

1.48 x 10~ 2 ° 

6 

0000 

0101 

0.0 

667. 379 

194.0 

8. 93 X 10' 18 

7 

0101 

0202 

667.379 

667.750 

15.0 

5.03 x 10" 19 

8 

0202 

0303(ll 3 0,03 3 0) n 

1335. 129 

668.180 

0.85 

2. 1 x 10 _2 ° 

9 

{10°0 > 02°0) I 1000 

1101(11 1 0,03 1 0) I 

1388.187 

688.678 

0.3 

7.0 x 10" 21 

10 

0101 

1000 (10°0, 02°0)j 

667.379 

720.808 

5.0 

1.7 x 10' 19 

11 

(llVoS 1 ®^ 0301 

1200(12°0, 04°0) I 

1932.470 

738. 643 

-u 

0.00832 

1.53 x 10 -22 

12 

0202 

1 1 

1101(11 0,03 0^ 

1335. 139 

741. 736 

0. 144 b 

3.52 x 10 -21 

13 

3 3 

(11 0,03 0) n 0303 

1202(12 2 0, 04 2 0) r 

2003.309 

757. 426 

0. 009 b 

2.0 X 10“ 22 

14 

(10°0,02 0 0) n 0200 

UOlfllVoS^)) 

1285,412 

791. 453 

0.022 

5. 5 X 10“ 22 

15 

(irV 03 * 0^0301 

1202(12 2 0,04 2 0) 1 

1932. 470 

828.265 

0. 0003 9 b 

-24 

7. 2 x 10 

16 

(10°0,02°0) I 1000 

0010 

1388. 187 

960. 955 

0.0184 b 

_oo 

4.38 x 10 

17 

(10°0, 02°0) n 0200 

0010 

1285.412 

1063.730 

0. 244 b 

6. 10 X 10' 22 


*The superscript b indicates the band intensity is from W. S. BENEDICT* < 20 ) the remaining intensities are from DRAYSON^ 21 ^ . 
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Figure 1 . Schematic of radiation components contributing to emergent 
radiation field at top of the atmosphere. 




(c) 


Figure 2. Schematic of variable wave number mesh for 0.1 cm" 1 inter- 
val with two sub-intervals in di . a) no molecular lines, 4 mesh points, b) 
one molecular line, 6 sub-intervals, 24 mesh points, c) two molecular 
lines, 11 sub-intervals, 44 mesh points. 
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Figure4. Water vapor continuum absorption coefficient for the 400-1400 cnr 1 region. The k 2 component 
is normalized to the self-broadened Modified Van Vleck-Weisskopf line shape for the rotational H 2 0 lines at 
400 cm' 1 and the self -broadened Lorentz line shape for the lines of the 1595 cm -1 H 2 0 band at 1400 cm'1 . 
The component is normalized in a similar manner to the foreign broadened water vapor components at 
400 and at 1400 cm'' . Only the ^ and k 2 contributions are included in the 400-1400 cm -1 region. 
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Figure 5. Homogenous path comparison for 1595 cm' 1 H 2 O band at room temperature, a) P { 

1 40 mb, P = 1 00 mb, p H 2 0 = 1 0 mb, u H 2 o = -0352 pr cm; b) P e = 1 07 3 mb, P = 985 mb, p H 2 o 
22 mb, and uh 2 o = -077 pr. cm. 
















